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SINGLE STAGE EXPERIMENTAL EVALUATION OF 
HIGH MACH NUMBER COMPRESSOR ROTOR BLADING 


PART I - DESIGN OF ROTOR BLADING 

uy 

D.R. Seyler and L.H. Smith, Jr. 


ABSTRACT 


A set of four axial-flow medium-aspect-ratio compressor rotor blade rows was 
designed to provide a vehicle for experimental evaluation of the use of blade 
camberline shape to minimize blade element losses in the high-transonic Mach 
number region. One blade row was designed for a tip diffusion factor of 0.35 
with a tip blade element ratio of supersonic to total camber of 0. The other 
three blade rows were designed for a tip diffusion factor of 0.45 with tip 
ratios of supersonic to total camber varying from 0 to 0.65, the value requir- 
ed for a double-circular arc blade element. 

SUMMARY 

A set of four axial-flow medium-aspect -rat io compressor rotor blade rows was 
designed to provide a vehicle for experimental evaluation of the use of blade 
camberline shape to minimize blade element losses in the high transonic Mach 
number region. One rotor, designated Rotor IB, was designed to produce a tip 
diffusion factor of 0.35 and the other three, designated Rotors 2B, 2D and 
2E, were designed to produce a tip diffusion factor of 0.45. Rotor IB de- 
livers a design total-pressure ratio of 1.60 and the Rotors 2 deliver a design 

total-pressure ratio of 1.76. The inlet hub-tip radius ratio is 0.50 for all 

four rotors and all were designed to pass a corrected weight flow of 29.66 
lb/sec-sq ft frontal area at a rotor tip speed of 1400 ft/sec. Radially con- 
stant total pressure but radially varying losses were used in the construction 
of velocity diagrams. The loss selections resulted in weight -flow-averaged 
rotor adiabatic efficiencies of 0.858 for Rotor IB and 0.837 for Rotors 2. 

The Rotors 2 were designed for three levels of the ratio of the camber in the 
front or supersonic portion of the blade element to the total camber. The 

smallest camber ratio was repeated in the design of Rotor IB. The aspect 

ratio for all is about 2.5; part -span shrouds are employed to insure satisfac- 
tory aeromechanical performance. The tip diameter of the compressor is 36.5”. 
For this diameter the design weight flow is 215.49 lb/sec. The compressor 
vehicle was designed to be compatible with the General Electric Lynn Component 
Test Operation House Compressor Test Stand. 

Design parameters were selected for inlet screening to permit the study of 
effects of inlet flow distortion. 
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INTRODUCTION 


Although early attempts to operate compressors with transonic Mach numbers 
were disappointing, the potential for substantial increases in stage pressure 
ratio with transonic compressors stimulated research in this area. In the 
early 1950* s the National Advisory Committee for Aeronautics embarked on an 
extensive program to develop the performance of transonic compressor stages 
culminating with the conclusion (ref, 1) that there is no Mach number barrier 
for axial-flow compressors, but rather a continuous spectrum of performance as 
design Mach number is increased into the supersonic range. 

This work and paralleling experiments in the aircraft gas turbine industry led 
to temporarily accepted design limits for compressor front end stages which 
acknowledged the observations that: 

1. Diffusion factor levels in excess of 0.4 in the rotor tip region are 
frequently accompanied by high losses. 

2. Blade elements with relative Mach numbers in excess of 1.2 have the 
potential for high losses when simple blade shapes such as the double- 
circular-arc are employed. 

3. Rotor tip solidities of less than 1.0 in conjunction with transonic 
relative Mach numbers accentuate the loss. 

Isolated demonstrations, such as that in reference 2, probed further with 
reasonable results, but sufficient data and a reasonable model for the design 
of front stage rotors with tip speeds in excess of 1200 feet per second do not 
appear in the literature. The purpose of the designs reported herein is to 
obtain blade element data from rotor blades which operate at high-transonic 
Mach numbers. 

Analyses of the origin of loss in transonic rotors, such as in reference 3, 
where the shock loss is considered separately from the profile or subsonic 
diffusion loss result in the observation that the shock loss can be a sizeable 
fraction of the total loss. Consideration of the relationship between the 
shock loss and supersonic expansion on the curved uncovered portion of the 
blade element led to the conclusion that the shape and direction of the suc- 
tion surface in the uncovered portion plays an important role in minimizing 
the loss of transonic rotor blade elements. 

It appears that, when high-transonic Mach numbers are employed, fundamentally 
different considerations apply to the selection of the shape of the uncovered 
forward portion of a blade element and the shape of the covered passage por- 
tion behind it. The uncovered portion, identified in sketch 1 as 1-m, may be 
further subdivided by the point, f, which represents the origin of the first 
Mach wave from the suction surface that is captured by the adjacent blade. 
Along the inlet portion, 1-f, expansion and compression waves, when they are 
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air flow 



sketch 1 


present, must be of such strength that the condition of no change in angular 
momentum of the flow along any blade-to-blade line forward of the leading 
edge is satisfied. A principal consideration in the selection of point, m, 
is the throat area requirement. This area is affected by starting conditions, 
streamtube contraction, boundary layer allowances and the change in choke 
flow area that occurs along a streamline because of changing radius in a ro- 
tating coordinate system. The shape of the covered portion, m-2 , of the blade 
depends upon subsonic diffuser considerations and deviation angle requirements. 

Finding the minimum blade element loss for high-transonic rotor blading seems 
to be most directly related to a search for a blade element shape which pro- 
vides the best balance between shock loss and subsonic diffusion loss. A new 
method for developing blade elements, wherein the camberline consists of two 
circular arcs which are mutually tangent at the point where they join, was 
employed. The front arc is identified as the supersonic arc, and the rear arc 
is identified as the subsonic arc. The term, camber ratio, refers to the 
ratio of the camber of the supersonic arc to the total camber. Blade elements 
developed in this way are called multiple-circular-arc elements and provide a 
reasonable means for varying blade shapes between the two extremes of double- 
circular-arc elements and multiple-circular-arc elements with the smallest 
supersonic camber consistent with choke-free operation. 
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For the 0.35 rotor-tip-diffusion-factor design the differences in camberline 
shape for the extremes approached normal manufacturing tolerances. Therefore 
only the blade shape approaching the choke limit was designed. For the 0.45 
rotor-tip-diffusion-factor design sufficient difference in camberline shape 
resulted so that the design of two extremes and an intermediate level of cam- 
ber ratio was warranted. Other design parameter selections are made so that 
the resulting blade geometry is typical aerodynamically and mechanically of 
future engine front stage rotor blade rows. 

The approach selected for these rotor designs is believed to offer the oppor- 
tunity for advancement in compressor technology by using: 

1. an axisymmetric procedure capable of calculating the flow field 
throughout the region of the blade row to accurately determine the 
flow direction at the blade edges and the distribution of meridional 
streamtube area through the blade row, and 

2. blade element shapes of the type described in appendix B, Multiple- 
and Double-Circular-Arc Properties, where the camberline and surfaces 
are circular arcs of different curvature in the uncovered or super- 
sonic portion and the covered or subsonic portion of the blade ele- 
ment . 

This multiple-circular-arc-blade type is believed to be capable of good per- 
formance if the proportions are properly selected. 


SYMBOLS 

The following symbols are used in this report: 

A exponent in equation (B12), page 26 

2 

A flow area, in 

a distance along chord line to position where maximum perpendicular dis- 

placement between camberline and chord line occurs, in 

enthalpy-equivalent stat ic-pressure-rise coefficient. 
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static -pres sure -rise coefficient , 
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blade chord length, in 

specific heat at constant pressure, 0.2399 Btu/lb-°R 

diffusion factor, equation (1), page 15 

base of natural system of logarithm 

2 

body force per unit mass, ft/sec 
frequency, cps 

acceleration due to gravity, 32.174 ft/sec 

displacement along radial direction between axisymmetric stream 
surfaces, in 

incidence angle, difference between air angle and camber line angle 
at leading edge in cascade projection (see Rotor Blade Shape Design 
section), deg 

mechanical equivalent of heat, 778.161 ft-lb/Btu 
total effective-area coefficient 

effective- area coefficient due to blade blockage 

effective-area coefficient due to part-span shroud blockage 

effective-area coefficient due to wall boundary layer blockage 

length in cascade projection (see Rotor Blade Shape Design section) 
defined in appendix C, in 

blade chord length projected to axial direction, in 
Mach number 

deviation coefficient, equation (2), page 7, and figure 9, page 54 

rotational speed, rpm 

number of rotor blades 

total or stagnation pressure, psia 
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p static or stream pressure, psia 

R gas constant, 53.342 ft-lb/lb-°R 

r radius, in 

r mean radius, average of streamline leading-edge and trailing -edge 

radii, in 

radius of curvature of projection of streamline on meridional plane, 
in 

S entropy, Btu/lb-°R 

s blade spacing, in 

T total or stagnation temperature, °R 

t static or stream temperature, °R 

t blade thickness, in 

t blade edge thickness, in; sketch 6 in appendix B 

e 

t blade maximum thickness, in 

m 

U rotor speed, ft/sec 

u circumferential displacement defined in equation (B19), page 30, in 

V air velocity, ft/sec 

w weight flow, lb/sec 

X empirical deviation angle adjustment, deg 

z displacement along compressor axis, in 

3 air angle, angle whose tangent is the ratio of tangential to axial 

velocity, deg 

Y ratio of specific heats 

y° blade-chord angle, angle in cascade projection (see Rotor Blade Shape 

Design section) between blade chord and axial direction, deg 

6° deviation angle, difference between air angle and camber line angle 

at trailing edge in cascade projection (see Rotor Blade Shape Design 
section), deg 
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£ meridional angle, angle between tangent to streamline projected on 

meridional plane and axial direction, deg 

0 angular displacement about compressor axis, deg 

n ad adiabatic efficiency 

1 unit vector in direction of intersection of ax i symmetric stream 

surface and blade mean surface 

—y 

< angle between cylindrical projection of i and axial direction, deg 

A lean angle, angle in z = constant plane between tangent to blade mean 

surface and radial direction, deg 

2 4 

P static or stream density, lb-sec /ft 

a solidity, ratio of chord to spacing 

$ camber angle, difference between angles in cascade projection (see 

Rotor Blade Shape Design section) of tangents to camberline at extremes 
of camberline arc, deg 

^ stream function; ip. = 0, ip =1 

h c 

w total-pressure-loss coefficient 

local loss-rate coefficient, equation (B3), page 24 
Subscripts: 

a point on camber line where maximum camber line rise occurs 

c tip or casing at any station 

cyl in cylindrical section of blade 

d design 

e equivalent two-dimensional cascade 

h hub at any station 

i arbitrary axial station 

i-1 preceding arbitrary axial station 

id ideal 
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m meridional direction 

m point on camberline where maximum thickness occurs 

r radial direction 

s suction surface 

sb subsonic or rear portion of blade element 

ss supersonic or front portion of blade element 

t tip at station 1.0 

t total when referring to blade element 

thr location in passage between adjacent blades where area bounded by 

adjacent blades and adjacent axisymmetric stream surfaces is minimum; 
appendix C 

z axial direction 

0 tangential direction 

1 leading edge 

1.0 z = constant plane indicating station at inlet to blade row 

1.5 z = constant plane indicating station at exit from blade row 

2 trailing edge 
Superscripts: 

* critical flow condition 

1 relative to rotor 


ROTOR DESIGN 
Overall Design Features 

A set of four medium-aspect ratio high-transonic rotors was designed as a 
vehicle for experimental evaluation of the use of blade camber line shape to 
minimize blade element losses. All four rotors have some overall character- 
istics in common as listed below: 

1. Rotor tip speed, 1400 ft/sec. 

2. Inlet hub-tip radius ratio, 0.50. 
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3. Rotor inlet tip diameter, 36,5 in. 

4. Corrected weight flow, 215.49 lb/sec. 

5. Corrected weight flow per unit frontal area, 29.66 lb/sec-sq. ft. 

6. Rotor tip solidity, 1.3 with chord radially constant. 

7. Number of rotor blades, 44. 

Two levels of loading result in differences in some of the overall character- 
istics as listed below: 



Rotor IB 

Rotors 2 

Rotor tip diffusion factor 

0.35 

0.45 

Total-pressure ratio, radially constant 

1.60 

1.76 

Rotor blade aspect ratio 

2.5 

2.4 

Rotor tip axial velocity ratio 

.91 

.91 

Rotor tip relative Mach number 

1.43 

1.45 

Tip blade element camber ratio 

1 level 

3 levels 

Weight -flow averaged rotor adiabatic efficiency 

0.858 

0.837 


The Rotor IB tip diffusion factor of 0.35 is representative of common practice 
for stages with an inlet radius ratio of 0.5; the tip Mach number, 1.43, is 
greater than is commonly used. This combination provides a link with past 
experience and other current designs where similar diffusion factor levels 
have been used in conjunction with lower Mach number levels. The selection of 
the Rotors 2 tip diffusion factor, 0.45, was made with the thought that the 
diffusion losses associated with a loading level which is somewhat higher than 
common practice might be tolerable, if the shock losses are minimized. The 
three levels of the ratio of supersonic to total camber which* characterize the 
difference between the three Rotors 2 are intended to provide a range inside 
which the minimum blade element loss will occur. One level of camber ratio 
for Rotor 2 was achieved by selecting double-circular-arc blade elements; the 
resulting blade is identified as Rotor 2D. The other two levels of camber 
ratio for Rotors 2 were achieved by using multiple-circular-arc blade elements. 
Rotor IB also employs multiple-circular-arc blade elements. All rotors employ 
double-circular-arc blade elements in the hub region. Because Rotor 2D is a 
blade composed entirely of double-circular-arc blade elements, it provides a 
geometric link with past experience. The remaining items were selected as 
being typical of a fan or compressor front stage, the most likely application 
of a high-transonic rotor blade. The radial position of the part-span shroud 
(app. A) is sufficiently far removed from the high Mach number region so as 
not to interfere with the important blade element measurements. The rotor tip 
diameter was selected to be compatible with the General Electric Lynn Component 
Test Operation House Compressor Test Stand. The overall flowpath is shown in 
figure 1. 


12 


Method of Calculation for Rotor Design 


The overall calculation procedure used for this design is outlined in Appendix 
B, Compressor Axisymmetric Flow Determination (CAFD), where the principal 
equations are listed. 

A total of twenty-three calculation stations and ten equal-flow streamtubes 
were used to arrive at the design vector diagrams in the region of the rotor. 
Figure 2 shows the location of calculation stations in the neighborhood of 
the Rotors 2. The first calculation station, 0.1, was located approximately 
15" forward of the rotor and the last calculation station, 2.0, was located 
approximately 15 n aft of the rotor (the annulus height at rotor inlet station 
is 9.125"); the streamline slope and curvature boundary conditions were es- 
tablished at stations 19" forward and 23" aft of the rotor. Two of the cal- 
culation stations were located immediately upstream and downstream of the 
blade axial projection extremities (stations 1.0 and 1.5, respectively) and 
five additional calculation stations were spaced equally within these ex- 
tremes. The additional interior stations improve the overall accuracy of the 
solution by including in fair detail the effects of blade thickness blockage 
on the slopes and curvatures at the blade inlet and exit stations and thereby 
enhance the execution of the rotor design. 

In the overall CAFD calculation the change in angular momentum is distributed 
along axisymmetric stream surfaces in the axial space of the blade according 
to the first quarter cycle of a sine wave raised to the A power (eq. (B12), 
page 26). In the present design, A was given the value of 1.5. Based on past 
experience this selection provides a reasonably good representation of the ex- 
pected rotor work input distribution; it is not a critical selection because 
the distribution does not have an important influence on the air angles at 
blade inlet and exit. Distributions of the change in angular momentum along 
two streamlines from Rotors 2 calculations are shown in figure 3. The dashed 
line in the region of the leading edge of the blade gives indication of the 
technique which is used to find modified axial derivatives of angular momen- 
tum that are used in the radial-equilibrium equation (eq. (B6) ). Modifying 
the derivatives minimizes the abnormal influence of large changes from point- 
to-point in the calculation grid. 

Effective-Area Coefficient 

In application of the CAFD procedure to the present designs, no attempt was 
made to calculate the localized velocity variations that occur deep in the 
annulus wall boundary layers. Instead, the calculated free-stream flow dis- 
tributions were continued to the annulus boundaries. The weight flow used at 
each station was related to the actual weight flow by an effective-area co- 
efficient, K, , , that accounts for the displacement thickness of the wall 

bkw 

boundary layers. Suitable values for the effective-area coefficient were 
selected on the basis of past experience; an effective-area coefficient of 
0.98 was used at all stations forward of the rotor and 0.95 at all stations 
aft of the influence of the rotor and the part- span shroud. 
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At calculation stations where physical obstructions to the flow exist, such as 
stations internal to a blade row, a modifier must be included in the effective- 
area coefficient to reflect the reduction of flow area. This modifier is a 
function of radius since both the circumference and the tangential thickness 
of the blades are dependent upon radius. The presence of part-span shrouds in 
these designs was reflected at the internal blade calculation stations by dis- 
tributing the resulting flow area reduction uniformly across the annulus; 
this method was judged most appropriate since an unusual and undesirable kink 
would have resulted from use of concentrated part-span shroud blockage. The 
blade thickness blockage was distributed radially as it actually occurred, 
however. The combined effects of boundary layer, blade thickness, and part- 
span shroud blockage resulted in minimum effective-area coefficients for the 
internal blade calculation stations varying from 0.85 at the tip to 0.73 at 
the hub for Rotors 2 and virtually the same values for Rotor IB. Figure 4 
shows the contributions of the three blockages which for a typical case make 
up the. total effective-area coefficient. As is the case with the angular 
momentum derivatives, the dashed lines in the region of the leading and trail- 
ing edge are used to find modified axial derivatives of total effective area 
coefficient. In these calculations the shroud blockage effect was applied 
over a larger axial distance than the actual shroud length. 

Annulus Shape and Axial Velocities 

For given overall design requirements, usually specified by weight flow, pres- 
sure ratio, and blade speed, the annulus area is the primary factor in deter- 
mining the level of axial velocity. When selecting axial velocity levels the 
following factors must be considered: (1) for a given value of blade loading 

parameter, such as diffusion factor, a reduction in axial velocity across the 
rotor tends to reduce the pressure ratio capability of the rotor, (2) a large 
increase in axial velocity tends to choke the annulus or to pose diffusion 
problems for succeeding blade rows, and (3) the slopes and curvatures of the 
casing and hub surfaces can sometimes be used effectively to yield axial 
velocity profiles that tend to strengthen the traditionally weak blade-end 
regions. An additional factor for consideration in this design is that the 
rotor is not immediately followed by a stator. Under this circumstance, 
annulus choke downstream of the rotor becomes an important consideration. For 
these designs it was decided that the tip axial velocity ratio should not be 
less than 0.9. 

The hub and casing contours for both Rotor IB and Rotors 2 annuli were estab- 
lished on the basis of the above considerations. In the course of the design 
calculations it became clear that it was possible to assign the same casing 
contour to both Rotor IB and Rotors 2 and still maintain a reasonable axial 
velocity change near the casing for both. Accordingly as is seen in figure 1 
only the hub contour is different for the two rotor types. Figure 2 is an 
enlargement of the flowpath in the neighborhood of the rotor. Included in 
this figure, for reference purposes, are the streamlines obtained from the 
design calculations for Rotors 2. 
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Total-Pressure Ratio and Loss Coefficients 


The selection of diffusion factor, solidity, inlet relative velocity and, from 
the preceding section, axial velocity ratio permits calculation for the rotor 
tip of the change in angular momentum across this section of the rotor. The 
change in angular momentum and a suitable rotor total-pressure-loss coefficient 
determine the rotor total-pressure ratio. The total-pressure ratio was main- 
tained constant over the radial height of the blade. The equations used for 
these calculations are given in reference 4, except the following definition 
of diffusion factor was used: 


D 


1 



r V 
2 02 


- r V 
1 01 


2roV^ 


( 1 ) 


This form is preferred when there is a change in radius as the flow passes 
through a blade row. 


The profile loss portions of the rotor total-pressure-loss coefficients used 
in the CAFD calculations were determined from the loss correlation (ref. 4) 
which is the object of additional development at NASA Lewis Research Center. 

In reference 4 it is demonstrated that losses from rotor tests for all radial 
immersions except the immersion ten-percent from the tip correlate reasonably 
well on diffusion factor when the total-pressure-loss parameter is used. Data 
at the ten-percent immersion above a diffusion factor of 0.25 scatters over a 
wide range of total-pressure-loss parameter. The designer has the choice of 
selecting values within this broad band of variation. The minimum level of 
the total-pressure-loss parameter band for ten-percent immersion was used to 
guide the selection of tip profile losses for all rotors. Additive shock 
losses were calculated by the NASA method given in reference 3. The diffusion 
factors and Mach numbers used to calculate the total-pressure-loss coefficients 
were taken from immediately preceding CAFD rotor design calculation and were 
therefore slightly different from those calculated in the final design. Fig- 
ure 5 shows the radial distributions of rotor total-pressure-loss coefficient 
that were employed. For Rotors 2 one loss curve was used to simplify the de- 
sign calculations. As will be shown in a later section reasonable radial 
distributions of adiabatic efficiency and diffusion factor resulted from 
these assumptions. 

Because it was decided to wash out the blockage influence of the part-span 
shroud downstream of the blade in a manner typical of the dissipation of a 
wake, it is necessary to provide some indication of the associated mixing 
loss. At the first calculation station, 1.54, where the part-span shroud 
blockage no longer exists a total -pres sure loss was introduced. The loss was 
determined by use of the compressible dumping loss coefficient associated with 
the two percent change in the effective area of the shroud and the compressi- 
ble velocity head at calculation station 1.5. The loss was distributed uni- 
formly from hub to tip. Accordingly, the part-span shroud loss did not in- 
fluence the shape of the rotor blading. 
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Rotor Blade Shape Design 


Constant -chord rotors employing multiple-circular-arc blade elements in the tip 
region and double-circular-arc blade elements in the hub region were selected 
to fulfill the requirements of these designs. The camberline of a multiple- 
circular-arc element consist of two circular arcs which are mutually tangent 
at the point where they join. The method of development of the multiple- 
circular-arc elements used in these designs is described in appendix B, 
Multiple- and Double -Circular -Arc Properties. The front arc is identified as 
the supersonic arc and the rear arc is identified as the subsonic arc. Geo- 
metric description of the double-circular-arc elements follows the customary 
procedure of using one circular arc to represent the camberline and circular 
arcs of different radii, tangent to the maximum thickness and edge thickness 
circles, to form the surfaces. The center of the maximum thickness circle is 
located on the camberline mid-way between the leading and trailing edges. 
Distributions of the ratio of supersonic to total camber are shown in figure 
6. Multiple-circular-arc-type elements extend inward to the radial location 
where the radii of the supersonic and subsonic arcs are equal. This radial 
location very nearly corresponds to the location of the part-span shrouds. 
Double-circular-arc -type elements are employed in the remaining radial length. 
Rotor 2D is composed of double-circular-arc-type elements along its entire 
radial length. 

The design of the rotor elements was performed along axisymmetric stream sur- 
faces using the projection recommended in reference 5. This projection cuts 
the blades along the axisymmetric surfaces but views the cut sections along a 
radial line. The incidence and deviation angles are defined in this projec- 
tion (identified as the cascade projection) and the resulting blade angles are 
termed cascade angles. 

The blade designs were based on the vector diagrams at the edges. In the hub 
region the inlet and exit calculation stations (1.0 and 1.5 in figure 2) are 
reasonably close to the leading and trailing edges. In the tip region however 
the edges are removed from the calculation stations by as much as twenty per- 
cent of the axial distance between the inlet and exit stations. The importance 
of properly accounting for this displacement will be seen later in this sec- 
tion. Rotor IB design data is listed in table 1; Rotors 2 design data are 
listed in tables 2, 3 and 4. Although only one overall CAFD calculation was 
performed for Rotors 2, individual vector diagrams were established for each 
Rotor 2 type depending upon the axial position of each in the annulus. These 
differences are observed when comparing the data in tables 2, 3 and 4. 

Design incidence angles were generally selected on the basis of past experience 
but subject to the additional criterion that the blade throat be adequate to 
preclude the condition of choking. The rotor geometry was checked for choking 
according to the method described in appendix C, Throat Area Parameter. In 
regions of elevated design inlet relative Mach number, it was intended that 
the design incidence angle be calculated according to the condition that the 


16 


suction surface angle at the leading edge be parallel to the upstream flow 
direction. When Rotor 2B was set according to this criterion and checked for 
choking, it was found that in the tip region less than 1% margin in throat 
area existed over the critical area associated with one normal shock at the 
upstream Mach number. In view of the uncalculated real flow effects this 
margin was judged to be insufficient to pass the design flow and accordingly 
the design incidence was increased. When the calculated design incidence for 
Rotor 2D was applied to Rotor 2B the margin shown in figure 7(a) resulted, and 
accordingly it was decided to use the Rotor 2D design incidence distribution 
to set the leading edge direction in the tip and pitch regions for all Rotors 
2. Similar results were encountered when Rotor IB tip blade elements were set 
to the desired incidence criterion and the design incidence for Rotor IB was 
likewise increased to that level which would have been used if a rotor had 
been designed with double-circular-arc sections to the loading level of Rotor 
IB. Figure 7(b) shows the final choke check results for Rotor IB. 

In the hub region, where the design inlet relative Mach number is near sonic, 
blade element losses are almost entirely related to diffusion requirements, 
and setting the suction surface at the leading edge parallel to the upstream 
flow direction is not a meaningful criterion. Design incidence angles were 
lowered to avoid excessive leading edge loading but limited by the additional 
consideration that the throat area be large enough to pass the design flow. 

The throat area ratio margin for the hub elements was selected to be about the 
same as the smallest margin for the tip elements. 


The tabulations of blade 
ences in direction, k 1 

* c I 


setting data (tables 1 through 4) list the differ- 
- , between the camber line and the suction surface 


at the leading edges. The incidence angles employed in these designs are 
shown in figure 8. 


The deviation angles were obtained by applying an extension of Carter’s Rule 
(fig. 9). The two curves in figure 9 for the positions of maximum camberline 
rise at forty and fifty percent of chord length are extracted from Carter’s 
correlation of cascade data in reference 6. The three other curves were con- 
structed by laying off at any blade-chord angle three increments in deviation 
coefficient equal to the increment between the forty percent and fifty per- 
cent curves. For each blade element an equivalent two-dimensional cascade 
was defined as is discussed in appendix B, Blade Setting. The deviation co- 
efficient, m, was obtained from figure 9 using the blade-chord angle and posi- 
tion of maximum camberline rise of the actual cascade. The deviation angle 
was calculated from the camber of the equivalent two-dimensional cascade. 
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The empirical adjustment, X, which is shown with the resulting deviation 
angles in figure 10 is derived from experience with aerodynamic design and 
performance synthesis for this general type of rotor. 

Some pertinent design data which summarize the results of the final overall 
CAFD calculations appear in figures 11 through 15 for Rotor IB and in figures 
16 through 20 for Rotors 2. Solid lines indicate edge conditions and broken 
lines indicate station conditions. Examination of the graphs gives indica- 
tion of the importance of properly accounting for the displacement between 
the stations and the edges. As an illustration, for Rotor IB the relative 
inlet air angle at the blade edge is up to two degrees different from that at 
calculation station 1.0 (fig. 13). The difference in design data (tables 2, 

3 and 4) for Rotors 2 are small enough so that with the scales selected 
they cannot always be clearly indicated in the graphs. Some indication of 
the design requirements for a stationary blade row to follow close onto the 
designed rotors is given in the graphs of absolute air conditions at the rotor 
exit calculation station (figs. 21 and 22). 

A layout on a cascade projection of a typical blade element from Rotor 2B 
is shown in figure 23. A portion of the pressure surface of the next blade 
in the counter-rotor-wise direction is included; this surface at any axial 
location is constructed with the correct circumferential displacement and its 
shape is therefore not identical with that of the first blade. The construc- 
tion of the throat width and upstream capture width which is performed to ar- 
rive at the Throat Area Parameter is included. Although this method of throat 
area determination is approximate, studies have shown that it is quite accurate 
except when the meridional angle, e, is very large and the number of blades is 
very small. 

The radial distributions of adiabatic efficiency (fig. 24) are obtained direct- 
ly from the overall CAFD calculations. But in keeping with the planned method 
for reduction of aerodynamic test data, where measured quantities will be cor- 
rected from the measurement planes to the blade edges, radial distributions of 
other performance indicators, diffusion factor, axial velocity ratio, relative 
turning angle, static-pressure-rise-coefficient and enthalpy-equivalent static- 
pressure-rise-coefficient (figures 25 through 29), were calculated from the 
design vector diagrams at the edges. 

The multiple-circular-arc and double-circular-arc blade elements were speci- 
fied in the cascade projection and the conversion from this projection to 
cylindrical sections was accomplished by means of the calculation routine 
which is discussed in appendix B, Special Airfoil Generator. Figure 30 is 
included to reveal the difference between cylindrical and cascade blade-edge 
angles for a typical case. In the hub region, differences as large as 3.5 
degrees occur, and thereby illustrate the importance of using the cascade 
projection. 
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For the purpose of manufacturing the blading, the developed cylindrical sec- 
tions are interpreted as plane sections normal to the stacking axis, since 
the error involved is negligible. 


TEST VEHICLE 
Inlet Ducting 

The conditions of zero inlet swirl and moderately thin casing boundary layer 
require incorporation of a flow straightener in the vehicle inlet and pro- 
vision for a large area contraction into the test vehicle; see figure 1. The 
flow straightener has honeycomb type cells which are approximately equivalent 
to a 7/8 m diameter; the length of the cells is 8". An incompressible axisym- 
metric flux plot of the inlet contraction region yields the casing velocity 
distribution shown in figure 31; the presence of the four struts was ignored 
for purposes of this calculation. The lack of adverse velocity gradients on 
both walls suggests conditions conducive to healthy boundary layers. 

Outlet Guide Vanes 

The same set of outlet guide vanes (figure 1) which was designed for the Con- 
tract NAS3-7618 vehicle (ref. 7) will be used for this vehicle. Since the 
performance of the rotor is of primary interest, the purpose of the outlet 
guide vane system is to remove swirl and decrease Mach number so that the ex- 
haust system pressure drop will be minimized when the discharge Mach number 
is high. The outlet guide vanes were designed to the calculations which were 
performed as a part of the Contract NAS3-7618 design. Check calculations 
were performed, and it was determined the outlet guide vane design performed 
under Contract NAS3-7618 would fulfill the needs of the test vehicle for 
Rotor IB and Rotors 2. Since the outlet guide vanes have variable stagger, 
any problems associated with their high incidence operation on the high- 
pressure part of the characteristic can be readily eliminated. 

The outlet guide vane sections employ a modified NACA 65-series thickness 
distribution and a circular-arc camberline. The hub solidity is 1.4, and 
the chord is constant radially. The aspect ratio is approximately 2. 

Distortion Screens 

Radial and circumferential distortion screens were designed to provide a 
total-pressure loss of 15% at the rotor inlet instrumentation station. The 
circumferential screen will cover a 90° arc of the inlet annulus while the 
radial screen will cover the outer 40% of the inlet annulus. The solidity of 
the screens for both types of distortion was determined to be 0.54. This 
value of solidity was arrived at utilizing test data of reference 8 which ap- 
proximates the desired design, and a correlation of screen total-pressure - 
loss coefficient with screen solidity and screen inlet Mach number (ref. 9). 
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A screen with approximately 0.016 inch wire diameter and 20 mesh will provide 
the required 0.54 solidity. These distortion screens will be supported by a 
screen covering the entire annulus whose wire diameter is about 0,092 inches 
and which has a 3/4 inch mesh. 


APPENDIX A 

ROTOR BLADE AEROMECHANICAL CONSIDERATIONS 

by 

P. Chifos and C.E. Danforth 


Summary 

The selected configuration for these rotors, titanium blading of aspect ratio 
of about 2.5 with part-span shroud at about 60-percent blade height, gives the 
most reasonable expectation of fulfilling the aeromechanical design objectives. 
Designs without part-span shrouds would have required aspect ratios as low as 
1.5. Although the reduced velocities* for the selected designs are in the 
range of reduced velocities for current blading operating successfully in 
General Electric engines, it is possible that self-excited blade vibration may 
be encountered in the 80-percent design speed region before stall is attained 
while the test vehicles are operated with circumferential inlet distortion. 

Resonant Vibration 

Two aeromechanical considerations governed the selection of the rotor blade 
designs : 

1. The avoidance of two per-rev resonance with the first flexural 
frequency of the blades at high speed. 

2. The attempt to guard against self-excitation as the incidence angle 
increases when stall is approached. 

While it is usually desirable to avoid integral order per-rev resonance in the 
low per-rev range and occasionally even higher order per-rev resonance, the 
avoidance of two per-rev resonance with the first flexural mode is generally 
mandatory for engine designs and for research vehicles which must be operated 
with circumferential distortion. Inlet systems for aircraft engines frequently 


* (Velocity relative to the blade/semi -chord times torsional frequency), dimen- 
sionless, varies approximately inversely with aspect ratio. 
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deliver one per-rev and two per-rev distortions, either of which can excite 
the first flexural vibration modes. The avoidance of two per-rev resonance 
may be achieved in either of two ways: 

1. Design for two per-rev flexural resonance at low speeds where ex- 
citation energy levels are correspondingly low, obtaining a blade 
with natural frequency at high speed between one- and two-per-rev. 
There are such designs which are currently successful. Aspect 
ratios are moderately high. 

2. Design to have a first flexural frequency above the two per-rev with 
adequate margin at high speeds, a result achieved by use of low as- 
pect ratios or of moderate-to-high aspect ratios with shrouding of an 
appropriate type. 

There are compromises in each approach. The low-speed resonance approach with 
correspondingly appreciable aspect ratios is workable only if the aerodynamic 
environment has been especially tailored to avoid those conditions conducive 
to instability in the comparatively high reduced velocity levels that natural- 
ly follow. The second approach of exceeding two per-rev frequencies demands 
a comparatively low aspect ratio (flexural frequency going up approximately 
as the chord) and a weight/size penalty. 

But the low aspect ratio achieving the flexural frequency margin over two 
per-rev carries with it a reduction in the reduced velocity which increases 
the design’s incidence migration range that is free from instability. How- 
ever, for the present aerodynamic designs this second approach would require 
an aspect ratio of about 1.5 with a conventional dovetail and no shrouding. 
Designs to the originally desired aspect ratio of three (3) imply the need 
for shrouding. 

After it became clear that part-span shrouds should be considered, parametric 
studies were made to assess the desirable spanwise location for the shroud. 
Figure 32 illustrates the choices available to avoid integral order resonance 
with the first flexural frequency. Resonance with two per-rev excitation is 
not a problem for any reasonable shroud location. The proposed Rotor IB de- 
sign yields a first flexural frequency at about 3.6 per-rev at design speed 
(1400 ft/sec tip speed, 8791 rpm) . Campbell diagrams are shown in Figures 33 
and 34 for Rotors IB and 2B, respectively. 

Self -Excitation 

A detailed examination of torsional frequency and corresponding reduced velo- 
city variation was made for aspect ratios 3.0 and 2.5. Primary emphasis was 
placed on Rotor IB which, because of its lower root camber and higher aspect 
ratio, has the lowest torsional frequencies. Whereas at 80-percent design 
speed, where the incidence margin appears to be the smallest, the aspect ratio 
3.0 design yielded a 1° incidence margin between clean inlet stall and self- 
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excitation, the aspect ratio 2.5 design results in a 3° incidence margin. Be- 
cause of the necessity to operate with inlet distortion and the unavoidable 
blade-to-blade manufacturing differences, the aspect ratio 3.0 design was con- 
sidered too risky. The aspect ratio 2.5 design was chosen. 

Steady Stress and Material Selection 

Calculations listed in Table 5 point clearly to the use of the titanium alloy, 
6A1-4V-Ti, which has the greatest service experience, on the basis of the fol- 
lowing considerations: 

1. Part-span shroud restraint and the resulting root warping stress to- 
gether with the stress components due to centrifugal and air loads 
combine to produce an untilted root edge stress of approximately 
100,000 psi for steel, practically the yield point for typical blad- 
ing material such as 403 stainless. This is an unreasonable choice in 
the event that any unanticipated vibration or malfunction be encoun- 
tered. Calculations indicate an untilted root edge stress of 69,000 
psi for titanium, which has a yield strength of 111,000 psi. This 
stress reduces to 13,500 psi for Rotor IB with tilting, as is shown 

in Table 5 . 

2. High tip speed and the related large blade pull makes the blade dove- 
tail design difficult with steel but relatively straightforward with 
titanium which has about the same strength but only 567o of the dead 
load . 

3. Elastic untwist and outer panel (the portion outside of the part-span 
shroud) stall stress without resonance appear to be acceptable. 


APPENDIX B 

DIGITAL COMPUTER PROGRAMS USED IN DESIGN 

Compressor Axisymmetric Flow Determination 

(CAFD) 

The name, Compressor Axisymmetric Flow Determination, and the abbreviation 
for it, CAFD, refer to a digital computer program which describes the axisym- 
metric flow field in the region of turbomachinery blading. 

The overall calculation scheme may be summarized as follows: The flow is 

considered at a number of axial stations (z = constant planes), and the radial- 
equilibrium equation, energy equation, and continuity condition are employed 
at each of them to determine the distribution of flow properties from hub to 
casing. It is necessary, however, that these distributions obtained separately 
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at each station be consistent from station to station and that the radial ac- 
celeration which a fluid particle undergoes as it passes from station to sta- 
tion be accounted for in the radial-equilibrium equation. This is done by 
assuming that the shape of a meridional streamline is adequately represented 
by a spline constrained at each axial station consistent with the continuity 
condition and at upstream and downstream boundary stations by a selected 
orientation and shape. The radial acceleration is expressed in terms of the 
meridional streamline slope and curvature. 

An iterative method of solution is implied. Meridional streamline shapes are 
assumed based on results from the previous iteration and flow distributions at 
each axial station are found. These imply new meridional streamline shapes, 
and this iterative calculation is allowed to continue until the changes in 
streamline shapes between two successive calculations are small enough to 
satisfy the user. 

An interesting feature of the computer program is the inclusion of terms in 
the radial-equilibrium equation which represent the blade action by a distri- 
buted body force field and the blade thickness by distributed blockage. The 
equation, which appears in reference 10, wherein physical interpretation of 
the individual terms is rendered, is restated here: 
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The heat addition term has been omitted since it does not appear in compressor 
design: in the absence of transients the blade surface and the fluid adjacent 

to it are at the same temperature. 

It is suggested (ref. 10) that a term be added to represent the heating effect 
of internal fluid friction, which encompasses all loss sources including 
shock waves : 
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The loss-rate coefficient, 


, is defined as follows: 
z 
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and is distinguished from the loss-coefficient, oj’ , in that the loss-rate 
coefficient is concerned with the manner in which the entropy increases in- 
side of the blade row whereas the loss coefficient merely acknowledges that 
the entropy has increased across the blade row. 

The radial and tangential components of the blade force may be expressed in 
terms of the axial change in angular momentum by means of the following re- 
lationships : 
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and 
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r 0 


(B5) 


By using the fact that the second derivative of the meridional streamline 

shape can be related to the meridional streamline curvature, — , and trigo- 

m 

nometry to relate velocity components, equation (Bl) may be rewritten as 
follows: 
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The particular grouping of terms on the right-hand side of equation (B6) is 
intentional: the first braces include terms which are present anywhere in the 

flow field while the second braces include terms which occur only when calcu- 
lations inside the axial extent of a blade row are desired. The attention of 
the reader is again invited to reference 10 where a numerical example illus- 
trates the relative magnitude of each term in equation (B6). 


At each axial station where the rotating blade row is influential the enthalpy 
of the fluid is related to the fluid swirl velocity by the energy equation 
which may be found in reference 11: 
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where 
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The static condition of the fluid is described by use of the equation of state, 

1 

p = gpRt , (B9) 

144 


and the specific heats are assumed to be constant. Application of the 
continuity condition in the form which may also be found in reference 11, 
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in combination with the previous four equations yields a complete description 
of the flow at each axial station. 


Because a form of the radial-equilibrium equation employing static pressure 
was used, an entropy gradient term does not appear explicitly. The entropy 
variations that result from radially non-constant losses are correctly ac- 
counted for, however, through the use of 
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which results from the definition of entropy and the assumption of a perfect 
gas. 
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With the flow field completely described, loading parameters, such as diffusion 
factor and static-pressure-rise coefficient, and forces on the blading and 
walls are calculated from station data using geometric characteristics of the 
system where required. 


Intra-Blade Energy Distribution 
(IBED) 


This digital computer program provides the stagnation pressure and temperature 
and the derivative of the product of radius and tangential velocity with re- 
spect to the axial direction along a streamline at desired axial locations in 
the region of a rotor blade row. The upstream and downstream stagnation pres- 
sure and temperature and the axial locations of the blade edges are specified. 
The axial variation of energy addition is calculated according to the follow- 
ing equation: 
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which allows analytical evaluation of the axial derivative of energy addition, 
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The stagnation temperature is determined from the energy equation, (B7), and 
the stagnation pressure is found by assuming a constant polytropic efficiency 
at the value associated with the downstream pressure and temperature. The 
value of A is selected from past experience to give a representative rate of 
energy addition consistent with expected blade surface velocity distributions. 

Blade Setting 
(BLAST) 

In this blade setting procedure, after the ax i symmetric flow field is estab- 
lished and the incidence and blade solidity are selected, the camberline 
orientation is calculated using Carter’s correlation (ref. 6). When changes 
in radius and axial velocity are present along a streamline across a blade 
row, the deviation angle is found by using the camber which is derived from 
an equivalent two-dimensional cascade of the same absolute circulation as is 
required along the streamline. 
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Sketch 2 shows a typical velocity diagram. The exit velocity diagram of the 
equivalent two-dimensional case appears in dashed lines. For the equivalent 
diagram to have the same circulation as the actual diagram, the equivalent ab- 
solute tangential velocity must be larger than the actual absolute tangential 
velocity by the ratio of the exit radius to the inlet radius. 



Addition of vectors in the tangential direction results in the following 
statement : 
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(B15) 


By substitution of the velocity diagram relationships, 

V = U - V tan s' , 

6 z 

and V z tan = V q » (B15a) 


and rearrangement, an expression relating the equivalent and actual relative 
fluid exit angles is obtained: 
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When stators are calculated the absolute fluid angle relationship is obtain- 
ed by setting the wheel speed to zero. 


Carter's formula, 
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may be expresed in terms of relative fluid angles and the incidence angle as 
follows : 
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The unguided portion of the fluid at the exit of the cascade is more closely 
related to the actual blade-chord angle than the blade-chord angle of the 
equivalent cascade and therefore the deviation coefficient, m, is determined 
by entering the functional relationship (shown graphically in figure 10) with 
the actual blade-chord angle. The trailing edge angle is calculated by ap- 
plication of the deviation angle from equation (B18) to the actual fluid exit 
angle. Iteration is required to arrive at the correct deviation coefficient. 
The iterative calculation is stopped when the blade-chord angle from the suc- 
cessive calculations differs by less than 0.01 degrees. 


Occasionally, past experience suggests the need for additional allowance in 
deviation angle, and an empirical adjustment, X, is provided in equation 
(B18) . The form used in this design is given in equation (2). 


Multiple- and Double-Circular-Arc Properties 

(MADCAP) 

For input to the Special Airfoil Generator computer program (described later 
in this Appendix) an accurate description of the blade camberline shape and 
thickness distribution is required. The MADCAP computer program provides for 
determination of the point at which two circular arcs join when multiple- 
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circular-arc blade elements are selected. In addition, the direction of the 
camberline and the circumferential displacement of the camberline are calcu- 
lated at preselected axial locations. After the trailing edge is located, the 
blade-chord angle is calculated. The angle between the camberline and the 
suction surface at the leading edge and the position along the chord line at 
which the maximum rise of the camberline occurs (in percent of chord) are also 
calculated. 

A section through a blade (sk. 3) is obtained by cutting the blade along a 
stream surface but viewing it in the radial direction as recommended in 





maximum 

thickness 


camber line 




sketch 3 



reference 5. This is called the cascade projection. The circumferential dis- 
placement of any point on the camberline from the leading edge in the cascade 
projection is defined as follows: 


u = 



(B19) 


In sketch 3 and other sketches in this section u is measured perpendicular to 
the compressor axis in the plane of the paper. 

The camberline of a multiple-circular- arc blade is made up of two circular arcs 
in this projection. The maximum thickness exists at the point where the two 
arcs join. The camber bisector of the forward or supersonic portion intersects 
the line of leading edges at a point midway between adjacent leading edges as 
shown. 


From knowledge of the inlet and exit air angles, the incidence angle and a 
tentative deviation angle, the leading edge direction, k| ,and the total camber 


are found. 


The value of camber ratio. 
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is specified and the total 


chord, c^ , and leading edge spacing, s^ , are known from preliminary design 


selections. The camber of the supersonic and subsonic portions are found sim- 
ply by multiplication and substraction. The chord of the supersonic portion 
may then be found from 
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By application of the law of cosines the subsonic chord may be found from 
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With the supersonic and subsonic chords calculated the radii of both arcs may 
be found from 
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When double-circular- arc blade elements are selected, specification of the 


value of camber ratio, 


ss 


, is not required; it can be determined from the 


*t 

observation that the radius of the camberline is the same in both the super- 
sonic and subsonic portions of the blade element, 


r ® r , 
ss sb 


2 sin — 


(B22b) 


and the use of equations (B20) and (B22): 
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The camber angles of the supersonic and subsonic portions may be found simply 
by multiplication and subtraction, and the supersonic and subsonic chords may 
be found from rearrangement of equations (B22) and (B22a) . 

The blade-chord angle of the blade element is the direction of the chord line, 
a straight line passing through the leading edge and trailing edge points, 
with respect to the axis of the compressor, and is found as follows; 
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Summation of the indicated lengths in the axial direction yeilds; 
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or rearranging: 
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The circumferential displacement from the leading edge may be determined in 
a similar manner: 
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To avoid numerical difficulty equation (B22) is substituted into equation 
(B27a) yielding the following preferred form: 
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and an alternate form is introduced for equation (B28) : 
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On the subsonic arc the following equations apply: 
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Again alternate equations are introduced to avoid numerical difficulty: 
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Determination of the chordwise position of the maximum camberline rise (per 
pendicular displacement from the chord line) depends upon the use of the 
observation that the direction of the camberline at the maximum camberline- 
rise point is the blade-chord direction. The axial and circumferential dis 
placements from the leading edge to the maximum camberline-rise point, a, a 
special cases of equations (B29) and (B30), where = y° . Therefore 



sketch 5 



Sketch 5 shows the case where the maximum camberline rise occurs toward the 
trailing edge with respect to the point, m, where the arcs join, or 

y° < (k ' - <t> ) . (B33) 

' 1 ss 

For the case where the maximum camberline rise occurs towards the leading 
edge with respect to the point where the arcs join, or 

Y° > (<; - <l> ss > , (B34) 

the axial and circumferential displacements may be calculated using equations 
(B27) and (B28) : 
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The displacement along the chord direction from the leading edge to the 
maximum camberline rise is obtained by addition of the two lengths indicated 
in sketch 5 and therefore 
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The difference in direction at the leading edge between the suction surface 
and the camberline is calculated following the method outlined in sketch 6. 

The edge thickness is considered to exist perpendicular to the camberline edge 
direction and at the edge; the suction and pressure surfaces are closed by a 
circular arc which passes through the camberline-edge point and is tangent 
to both surfaces. On the maximum thickness line the suction surface arc and 
the camberline arc have the same direction. The incidence angle required to 
set the suction surface at the leading edge parallel to the upstream flow 
direction is therefore just the difference in camber in the supersonic portion 
between the suction surface arc and the camber line arc: 
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the case in the tip region for Rotor 2D) ; the camber of the arc between the 
leading edge and the maximum thickness points can then be designated (f> ^ , and 

the chord of the arc, c , in both sketch 6 and equation (B40) . By observing 

m 

that for double-circular- arc elements the following relationships exist: 


*m 2 


and 


c - 
m 


2 cos 


m 


(B41) 


for the double-circular arc equation (B40) may be written 


i = 2 arctan 


Y t Y t 

c tan 7 — + t - t cos — 
t 4 m e 2 


(B42) 


c + t sin — 
t e 2 


Special Airfoil Generator 
(SAG) 


With this program the description of the blade as a series of cascade projec- 
tions for the several axisymmetric stream surfaces is converted to a descrip- 
tion on cylindrical surfaces from which the blade can be manufactured. The 
blade is stacked and lean angles, X , are calculated in the event additional 
axisymmetric flow field calculations are desired. 


At each grid point (intersection of meridional streamline and axial station) 
cylindrical blade angles are calculated from 

K ' C yl = arctan (tan K f - tan e tan X) , (B43) 

and curve-fit radially so that values may be found on the desired cylindrical 
surfaces . 


From knowledge of the axial locations of the blade edges, edge angles are cal- 
culated and an analytical description of the camberline angle versus z on the 
cylindrical surface is obtained. The effective-area coefficient due to the 
blade blockage and the edge thickness are curve-fit radially and values on the 
desired cylindrical surfaces are found. Coordinates for the camberline are 
obtained by integrating the camberline angle equation using additional axial 
locations for more complete description of the blade section. Surface coordi- 
nates are found by calculating the thickness distribution from 


t 



(1 - W 2,r cos K cvl 


Vt 


(B44) 


and applying it normal to the camberline. 
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After the centers-of-area of the sections are determined the sections are 
displaced axially and circumferentially to obtain the desired relationship 
of each section to a radial line. Lean angles are calculated at the desired 
axial planes by differentiation of the functional relationship between circum- 
ferential position and radius. The sign convention for lean angle, A , is 
that which yields a positive angle when a point moving along the blade moves 
opposite to the direction of rotor rotation as it moves outward along the 
blade. 


APPENDIX C 


THROAT AREA PARAMETER 


The knowledge of the details of the axi symmetric flow field inside a tran- 
sonic rotor blade row allows an approximate calculation which measures the 
capability of any rotor blade element to pass the flow for which it was de- 
signed. The result of this calculation is named the Throat Area Parameter; 
it is defined as follows: 


thr _ 


^bkthr^thr ^thr 


(Cl) 


* 

A i 


thr 


The radial stream tube height, h, is obtained from overall CAFD calculations. 
The capture width, , is just the product, s^cos 6 ? - The throat width, 


L is the distance between the pressure and suction surfaces of a cascade 
thr 

in the cascade projection at the location where the area found by the product, 
Lh ... 

— , is minimum. An example of the construction used to provide informat- 

A thr 

ion for Throat Area Parameter calculations appears in figure 23. Trial-and- 
error is required to arrive at the throat location. The pressure surface of 
the next blade in the counter-rotor-wise direction is constructed with the 
correct circumferential spacing and its shape is therefore not identical with 
that of the first blade. 


The effective area coefficeint, K , , , is introduced as a modifier on the 

bkthr 

throat stream tube height to provide for inclusion of any influence on the 
throat area which is not otherwise introduced. In these designs the blockage 
of the part- span shroud was distributed equally over the annulus height and 
accordingly resulted in no additional contraction of the axisymmetric stream- 
tube. The two percent blockage (fig. 4) of the part span shroud was used in 
all calculations (figs. 7(a) and 7(b)) involving equation (Cl). 
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The change in critical flow area between the leading edge and the throat re- 
sults from the change in relative stagnation conditions as the streamline 
radius changes. The ratio of critical areas may be expressed as follows: 

i±i 


thr 


P thr 

p ; 


2y 


id 


(C2) 


where the ideal relative stagnation-pressure ratio is given by a rearrange- 
ment of an equation found in reference 11, 



si + 


l^l) 

2 | 

' izl 

V' j 


l 2 . 

1 11 


1 1 + 



2 




Y-l 


(C3) 


The Throat Area Parameter is usually plotted against inlet relative Mach num- 
ber as in figures 7(a) and 7(b). On these coordinates, curves of throat area 
required for one-dimensional isentropic flow and for flow with the loss in 
relative stagnation pressure associated with one normal shock at the inlet 
relative Mach number may be constructed. A Throat Area Parameter in excess 
of the one-normal-shock curve is preferred because the procedure for checking 
the throat area does not account for the following real flow effects which 
tend to require more throat area: 


1. build-up of blade boundary layer, 

2. build-up of annulus wall boundary layers from the leading edge to 
the blade throat, 

3. shock losses in excess of those associated with one norma? shock at 
the inlet relative Mach number, and 

4. lack of uniform flow in the free stream at the blade throat. 
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Table 1. - Cascade Projection Data for Rotor lB Blade Setting 
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Table 2. - Cascade Projection Data for Rotor 2B Blade Setting. 
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Table 3. - Cascade Projection Data for Rotor 2E Blade Setting 
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Table 4. 

- Cascade Projection Data 

for Rotor 2D Blade 

Setting. 
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Table 5. - Aeromechanical Data for Rotor Blading 


Rotor IB 


Rotor 2B 


Design speed stresses (psi) 

Average root centrifugal stress 29, 600 

Total stress at root 

Convex side 59, 500 

Trailing edge 13,500 

Total stress at shroud 

Convex side 16, 900 

Trailing edge 11,000 


Maximum estimated stall vibratory 
stress at shroud (7 0 endurance) 

Convex side 8,400 (22) 

Trailing edge 7,600 (18) 


28, 600 


27,000 
27, 000 


15, 500 
12, 500 


12, 300 (31) 
5, 300 (13) 


Blade natural frequencies at design speed (cps) 


1st flexural 
2nd flexural 
1st torsional 

Tip untwist at design speed (deg ) 
D esign speed reduced velocity 
Material 


525 

575 

1910 

2150 

1205 

1340 

0.8 

0.7 

1.33 

1.19 

6A1-4V-Ti 

titanium alloy 
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FIGURE I. COMPRESSOR FLOWPATH 
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AXIAL DISTANCE ,Z, IN 


FIGURE 2. DESIGN FLOW FIELD FOR ROTORS 2 
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Station 0.98 1.0 1.017 1.05 1.083 1.51 

1.033 1.067 1.3 



Axial distance, z , in 

Figure 3. - Typical axial variation of angular momentum. 




Figure 4. - Typical axial variation of effective area coefficient. 
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Relative total-pressure-loss coefficient 
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Figure 5. - Radial variation of design relative total-pressure-loss coefficient 
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Figure 7(a). - Rotors 2 design choke check. 
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Relative Mach number, Mj 


Figure 7(b). - Rotor IB design choke check 
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Deviation coerncient, 



Blade-chord angle , Y ° , deg 


Figure 9. - Coefficients for design deviation angle rule. 
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Deviation angle 



Figure 10(b). - Radial variation of Rotors 2 design deviation angle 

and empirical deviation angle adjustment. 
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Figure 11. - Radial variation of Rotor IB design axial velocity. 
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Relative Mach number, 



Figure 12 


Radial variation of Rotor lB design relative Mach number 
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Figure 13 


Radial variation of Rotor IB design relative air angle 
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Figure 14. - Radial variation of Rotor IB design streamline meridional angle 
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Figure 15. - Radial variation of Rotor IB design streamfuction. 
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Axial velocity, V , ft/sec 
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Figure 16. - Radial variation of Rotors 2 design axial velocity, 
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Figure 17 


Radial variation of Rotors 2 design relative Mach number 
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Figure 18. - Radial variation of Rotors 2 design relative air angle. 
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Figure 19. - Radial variation of Rotors 2 design streamline meridional angle 
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Figure 20. - Radial variation of Rotors 2 design streamfuction, 
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Figure 25. - Radial variation of design diffusion factor. 
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Figure 29. - Radial variation of design enthalpy -equivalent 
static -pres sure -rise coefficient. 
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Figure 31. - Surface velocity distribution through inlet contraction region, 
(Incompressible axi symmetric solution) 
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IF = 1st. flexural frequency 
IT = 1st. torsional frequency 
2F = 2nd. flexural frequency 
* = design speed (8791 rpm) 


Figure 33. - Rotor IB Campbell diagram. 
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* = design speed (8791 rpm) 

Figure 34. - Rotor 2B Campbell diagram. 
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